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We apply new upper limits on neutrino fluxes and the diffuse extragalactic component of the GeV 
7— ray flux to various scenarios for ultra high energy cosmic rays and neutrinos. As a result we find 
that extra-galactic top-down sources can not contribute significantly to the observed flux of highest 
energy cosmic rays. The Z-burst mechanism where ultra-high energy neutrinos produce cosmic rays 
via interactions with relic neutrinos is practically ruled out if cosmological limits on neutrino mass 
and clustering apply. 



I. INTRODUCTION 

High energy neutrino astrophysics is currently very 
active, in particular experimentally [lj. Neutrino tele- 
scopes are reaching sensitivities comparable to theoreti- 
cal expectations for neutrino fluxes based on their con- 
nection to primary cosmic rays and secondary 7— rays 0- 
This is in particular the case for upper limits from 
AMANDA II in an energy range between ~ 10 14 eV 
and ~ 10 18 eV, and from the RICE experiment Q above 
10 16 eV. The former aims to detect neutrinos by looking 
for showers and/or tracks from charged leptons produced 
by charged current reactions of neutrinos in ice, whereas 
the latter is searching for radio pulses emitted by neu- 
trino induced showers in south polar ice. In addition, 
based on the non-observation of radio pulses from the 
Earth's surface expected from neutrinos above ~ 10 22 eV, 
the FORTE satellite has established upper limits on their 
fluxes in this hitherto unexplored territory p| . 

In the near future sensitivities will further improve 
by next generation versions of these techniques. For 
AMANDA this will be ICECUBE [(J at the South pole 
and possibly a comparable kilometer scale neutrino tele- 
scope in the Mediterranean , based on ANTARES 
NEMO 0, and NESTOR [<j. Improved limits from the 
radio technique may come from the Antarctic Impulsive 
Transient Antenna (ANITA) which is a planned long du- 
ration balloon mission to detect radio waves from showers 
induced by neutrinos in the antarctic ice [Io| . 

Next generation experiments for ultra-high energy cos- 
mic rays (UHECR) above ~ 10 19 eV will also have consid- 
erable sensitivity to neutrinos, typically from the near- 
horizontal air-showers that are produced by them 
These projects include the southern site of the Pierre 
Auger Observatory |12| , a combination of a charged par- 
ticle detector array with fluorescence telescopes for air 
showers produced by cosmic rays above ~ 10 19 eV, and 
the telescope array ^3|, which may serve as the opti- 
cal component of the northern Pierre Auger site. There 
are also plans for space based observatories such as 
EUSO [3 and OWL |l7| of even bigger acceptance. 



Finally, there are plans to construct telescopes to de- 
tect fluorescence and Cerenkov light from near-horizontal 
showers produced in mountain targets by neutrinos in 
the intermediate window of energies between ~ 10 15 eV 
and ~ 10 19 eV 0,H3- Acoustic detection of neutrino 
induced interactions is also being considered |2l| . 

In an earlier paper p3 | we reviewed fluxes in various 
scenarios in the context of constraints from current cos- 
mic ray data and upper limits on 7— ray and neutrino 
fluxes. Besides the improved neutrino flux limits from 
AMANDA II, RICE, and FORTE, a possibly lower ex- 
tragalactic contribution to the diffuse GeV 7— ray back- 
ground observed by the EGRET instrument on board 
the Compton 7— ray observatory piiL I24I |25| has been 
pointed out recently. An upper limit on the extragalactic 
diffuse 7— ray flux constrains the total amount of electro- 
magnetic (EM) energy injected above ~ 10 15 eV which 
cascades down to below the pair production thresh- 
old for photons on the cosmic microwave background 
(CMB) |2a.l27|. Since in any scenario involving pion pro- 
duction the EM energy fluence is comparable to the neu- 
trino energy fluence, a change in the constraint on EM en- 
ergy injection can also influence allowed neutrino fluxes. 
Furthermore, future 7— ray detectors such as GLAST |28| 
will test whether the diffuse extragalactic GeV 7— ray 
background is truly diffuse or partly consists of discrete 
sources that could not be resolved by EGRET. This could 
further improve the cascade limit. 

Motivated by these improved constraints and 
prospects, and by more detailed information available 
on neutrino sensitivities of future experiments, in the 
present paper we reconsider flux predictions in scenar- 
ios currently often discussed in the literature. As in 
Ref. l22l. we apply our recently combined propagation 
codes [2JJ, |3(j, |3ll I33] ■ Sect. II summarizes the numerical 
technique used in this paper. In Sect. Ill we discuss the 
cosmogenic neutrino flux, i.e. the flux of neutrinos pro- 
duced as secondaries of extragalactic cosmic rays during 
propagation, and its dependence on various UHECR 
source characteristics. In Sect. IV we review neutrino 
flux predictions in extragalactic top-down scenarios 
where UHECRs are produced in decays of super-massive 
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particles continuously released from topological defect 
relics from the early Universe. If UHECRs are new 
hadrons, they have to be produced as secondaries of 
accelerated protons which also gives rise to neutrinos. 
Their fluxes in these scenarios are reviewed in Sect. V. 
Sect. VI discusses primary neutrino fluxes required in 
scenarios where the cosmic rays observed at the highest 
energies are produced as secondaries from interactions 
with the relic cosmological neutrino background, often 
called Z-burst scenario. In Sect. VII we discuss neutrino 
fluxes from active galactic nuclei (AGN) models. Finally, 
in Sect. VIII we conclude. 



III. THE COSMOGENIC NEUTRINO FLUX 

The flux of "cosmogenic" neutrinos is created by decay- 
ing charged pions produced in interactions of primary nu- 
cleons of energy above ~ 5 x 10 19 eV with CMB photons, 
the Greisen-Zatsepin-Kuzmin (GZK) effect |37j. This 
flux depends on the production rate of the primary nu- 
cleons which we parameterize according to Eq. ft}. 



A. Dependence on diffuse photon background 



II. NUMERICAL TECHNIQUE 

Our simulations use an implicit transport code that 
evolve the spectra of nucleons, 7— rays, electrons, 
electron-, muon-, and tau-neutrinos, and their anti- 
particles along straight lines. Arbitrary injection spec- 
tra and redshift distributions can be specified for the 
sources and all relevant strong, electromagnetic, and 
weak interactions have been implemented. For details 
see Refs. [2^, H^, |3(], HH- Our results apply to aver- 
age large scale extragalactic magnetic fields of the order 
B < 1CP 11 G and a universal radio background between 
the minimal values consistent with observations |33l ] and 
modera tely low theoretical estimates from radio source 
counts j2J|. 

For the neutrinos we assume for simplicity that all 
three flavors are maximally mixed which for our purposes 
is an excellent approximation |35ll36j| and thus have equal 
fluxes. For each flavor we sum fluxes of particles and 
anti-particles. 

In the present investigation we parameterize power law 
injection spectra of either protons (for UHECR sources) 
or neutrinos (for Z-burst models) per co-moving volume 
in the following way: 

<j>(E,z) = f(l + z) m E- a e(E m ^-E) 

^min — Z ^ Z maK , (1) 

where / is the normalization that has to be fitted to 
the data. The free parameters are the spectral index 
a, the maximal energy £ mM , the minimal and maximal 
redshifts z m j n , z max , and the redshift evolution index m. 
The resulting neutrino spectra depend insignificantly on 
z m i n in the range < z m ; n < 0.1 where local effects could 
play a role, and thus we will set z m - ln = in the following. 

To obtain the maximal neutrino fluxes for a given set 
of values for all these parameters , we determine the max- 
imal normalization / in Eq. Q by demanding that both 
the accompanying nucleon and 7— ray fluxes are below 
the observed cosmic ray spectrum and the diffuse 7— ray 
background observed by EGRET, respectively. 
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FIG. 1: Dependence of the average cosmogenic neutrino flux 
per flavor (labeled "1/4") on the contribution of the accompa- 
nying photon flux (labeled "7" ) to the old (upper error bars on 
the left) and new (lower error bars on the left) diffuse EGRET 
flux estimate. The proton primary parameters in Eq. Q have 
been fixed to m = 3, z max = 2, and a = 1, whereas -E max was 
varied. The UHECR proton flux (labeled " p" ) is normalized 
to the data of the AGASA |3g| and HiRes |39t experiments. 
The numbers indicate the fraction of the new EGRET esti- 
mate of the extragalactic diffuse 7— ray flux contributed by 
the respective scenario, where the unlabeled thick curves cor- 
respond to 100%. 

We first consider the dependence of the cosmogenic 
neutrino flux on the contribution of the accompanying 
photons to the diffuse 7— ray flux in the lOOMeV - 100 
GeV region. As an example, we do this by fixing m = 3, 
Zmax = 2, and a = 1 in Eq. ft}, while varying the max- 
imal energy E m&K . The relatively hard proton spectrum 
in this scenario could be produced, for example, by ac- 
celeration in potential drops or reconnection 40]. Wo 
normalize the resulting UHECR flux to the observations 
above 10 19 eV and note that the discrepancy between the 
AGASA H3 and HiRes (H fluxes above ~ 10 20 eV has a 
negligible influence on the predicted cosmogenic neutrino 
flux. As Fig. n shows, a decrease of the diffuse photon 
flux results from a decrease of -Emax in this case. The 
old EGRET flux estimate corresponds to E max = 
3 x 10 22 eV, whereas the ~ 50% smaller new EGRET 
flux corresponds to E max = 2 x 10 22 eV. Most likely 
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only a fraction of the measured diffuse photon back- 
ground is connected to GZK neutrinos. In the present 
scenario 0.5 and 0.2 of the flux measured by EGRET 
corresponds to E max = 10 22 eV and E ms , x = 3 x 10 21 eV, 
respectively. Note that the UHECR proton flux is the 
same in all cases, except for the highest not yet observed 
energies where the flux can be affected by the distance 
to the nearest sources. 

The cosmogenic neutrino flux has recently been re- 
evaluated also in Ref. |4l| where the EGRET constraint 
has not been taken into account. The latter, however, 
eliminates a considerable part of the higher fluxes con- 
sidered there. 



B. Comparison with experimental limits and 
future sensitivities 
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FIG. 2: A scenario with maximal cosmogenic neutrino fluxes 
per flavor as obtained by tuning the parameters of the proton 
primaries in Eq. Q to 2: max = 2, -E max — 2 x 10 22 eV, m = 3, 
a = 1. Also shown are predicted and observed cosmic ray 
and 7— ray fluxes, the atmospheric neutrino flux |42|. as well 
as existin g u pper limits on the diffuse neutrino fluxes from 
MACRO gjj, AMANDA II 0L BAIKAL 0, AG AS A 0, 
the Fly's Eye |4?j | and RICE (4J experiments, and the limits 
obtained with the Goldstone radio telescope (GLUE) | 4q | and 
the FORTE satellite |f|, as indicated. The cosmic ray data 
are as in Fig. whereas only the new EGRET flux is shown 
to the left. 

The two major categories of experiments are based on 
detection in water, ice or underground, typically sensi- 
tive below ~ 10 16 eV, and on air shower detection, usu- 
ally sensitive at higher energies. Existing neutrino flux 
upper limits come from the underground MACRO ex- 
periment P3 at Gran Sasso, AMANDA II in the 
South Pole ice, and the Lake BAIKAL neutrino tele- 
scope 1441 in the first category, and the AGASA ground 
array |45l . the former fluorescence experiment Fly's 
Eye [4y|, the Radio Ice Cerenkov Experiment RICE H 
(there is also a limit from the HiRes experiment |49| 
which is between the RICE and AGASA limits), the 
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FIG. 3: The cosmogenic neutrino flux per flavor shown 
in Fig. in comparison with expected sensitivities of 
the currently being constructed Pierre Auger project to 
tau- neutrinos [Tl|. the planned projects telescope a rray 
(TA) |B(J, the fluorescence/Cerenkov detector NUTEL l2(fl. 
the space based EUSO pill , the water-based Baikal [44j| 
and ANTARES |2 (the NESTOR sensitivity for 1 tower 
would be similar to AMANDA-II and for 7 towers sim- 
ilar to ANTARES 0), the ice-based AMANDA-II 
and ICECUBE (similar to the intended Mediterranean 
km 3 project 0), and the radio detectors RICE (5^1 and 
ANITA [h|, as indicated. All sensitivities except for ANITA 
and RICE refer to one year running time. For comparison, the 
7— ray bound derived from the EGRET GeV 7— ray flux Q 
is also shown. 



Goldstone Lunar Ultra-high energy neutrino experiment 
GLUE 0, and the Fast On-orbit Recording of Tran- 
sient Events (FORTE) satellite H in the second cate- 
gory. As an example, an optimistic cosmogenic neutrino 
flux is compared with current neutrino flux upper lim- 
its in Fig. [3 Future experiments in the first cate gory 
include NT200+ at Lake Baikal 0, ANTARES |H2|, 
NESTOR in Greece 0, as well as a possible common 
km 3 scale detector in the Mediterranean 0, HI , and ICE- 
CUBE 6], the next-generation version of AMANDA at 
the South pole. The air shower based category includes 
the Pierre Auger project the telescope array [o0| . 
the fluorescence/ Cerenkov detector NUTEL [20|. and the 
space based EUSO [EJ and OWL [55| experiments. The 
EUSO sensitivity estimate used here is based on deeply 
penetrating atmospheric showers induced by electron or 
muon-neutrinos only |5lj| and may thus be considerably 
better if tau neutrinos, Cerenkov events, and Earth skim- 
ming events are taken into account |56|. for which there 
are no final estimates available yet. The same applies 
to the OWL project [5{|. The cosmogenic neutrino flux 
models shown in Fig. ^ are compared with future sensi- 
tivities in Fig. |31 

The fluxes shown in Figs. |3 and 0| are consider ably 
higher than the ones discussed in Refs. [oftl5&lofllr^loll |. 
and should be easily detectable by at least some of these 
future instruments, as demonstrated by the expected 
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FIG. 4: The cosmogenic neutrino flux average per flavor as 
shown in Fig. in comparison with differential sensitivities 
expected for 2006 for the Pierre Auger project [TH (for two 
years of data taking, assuming strong (upper curve) and no 
(lower curve) deep inelastic tau lepton scattering) and the 
ANITA project 10] for 10 (upper curve) and 30 (lower curve) 
days of observation. The corresponding number of events ex- 
pected are listed in Tab.[j] 



Experiment 


F„(1.0) 


F„(1.5) 


F„(0.5) 


F,(0.2) 


2 year Pierre Auger min 


4.5 


4.8 


3.2 


1.8 


2 year Pierre Auger max 


13.4 


14.5 


9.6 


5.4 


10 days ANITA 


5.0 


5.8 


2.9 


1.3 


30 days ANITA 


14.9 


17.5 


8.7 


3.8 



TABLE I: The number of tau neutrino (Pierre Auger) and 
electron neutrino (ANITA) events expected to be measured 
by 2006 for the four neutrino fluxes shown in Figs. [3] and 2] 
The Pierre Auger differential sensitivity shown assumes two 
years of data taking until 2006. The minimal (upper line in 
Fig.|IJ and maximal (lower line in Fig.|IJ sensitivity depends 
on the assumptions made for the strength of deep inelastic 
scattering of tau leptons. For ANITA the case for balloon 
flights of 10 and 30 days are shown. 

number of events listed in Tab. [H 

We note that the non-observation of GZK neutrinos in 
2006 will significantly restrict the accompanying photon 
contribution to the EGRET diffuse flux. 



IV. NEUTRINO FLUXES IN TOP-DOWN 
SCENARIOS 

Historically, top-down (TD) scenarios were proposed 
as an alternative to acceleration scenarios to explain the 
huge energies up to 3 x 10 20 eV observed in the cosmic 
ray spectrum [62]. In these top-down scenarios UHE- 
CRs are the decay products of some super-massive "X" 
particles of mass mx 3> 10 20 eV close to the grand uni- 
fied scale, and have energies all the way up to ~ mx- 
Thus, the massive X particles could be meta-stable relics 



of the early Universe with lifetimes of the order the cur- 
rent age of the Universe or could be released from topo- 
logical defects that were produced in the early Universe 
during symmetry-breaking phase transitions predicted by 
in Grand Unified Theories (GUTs). The X particles 
typically decay into leptons and quarks. The quarks 
hadronize, producing jets of hadrons which, together 
with the decay products of the unstable leptons, result in 
a large cascade of energetic photons, neutrinos and light 
leptons with a small fraction of protons and neutrons, 
some of which contribute to the observed UHECR flux. 
The resulting injection spectra have been calculated from 
QCD in various approximations, see Ref. [27] for a review 
and Ref. [631 ] for more recent work. In the present work 
we will use the QCD spectra discussed in Ref. and 
shown in Fig. 11 of Ref. [22j. For the purposes of the 
current work this is not expected to make a significant 
difference as compared to the more accurate fragmenta- 
tion spectra |65| . 
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FIG. 5: Flux predictions for a TD model characterized by 
p = 1, mx — 2 x 10 13 GeV. The contribution of photons and 
protons to the UHECR flux decreases with decreasing frac- 
tional contribution to the diffuse photon flux at EGRET en- 
ergies which is denoted in numbers. Even a TD contribution 
to the present estimate of the diffuse EGRET flux as high as 
100% (unlabeled thick curves) is only marginally consistent 
with the AGASA UHECR excess. The line key is the same 
as in Fig.0 

For dimensional reasons the spatially averaged X par- 
ticle injection rate can only depend on the mass scale mx 
and on cosmic time t in the combination 

h x {t) = K m p x t~ 4+P , (2) 

where k and p are dimensionless constants whose value 
depend on the specific top-down scenario |62j |. Extra- 
galactic topological defect sources usually predict p = 1, 
whereas decaying super-heavy dark matter (SHDM) j6|| 
\fil\ of lifetime much la rger than the age of the Universe 
corresponds to p = 2 [23. It has also been suggested 
that annihilation of SHDM particles instead of their de- 
cay might contribute to the observed UHECRs [6^ ■ 
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In the SHDM scenario the observable UHECR flux will 
be dominated by the decay or annihilation products of 
SHDM in the Galactic halo and thus by sources at dis- 
tances smaller than all relevant interaction lengths. Com- 
position and spectra will thus be directly given by the 
injection spectra which are dominated by photons. This 
is most likely inconsistent with upper limits on the ultra- 
high energy (UHE) photon fraction above 10 19 eV [69$ . 
However, requiring that this scenario explains only UHE- 
CRs above 4 x 10 19 eV allows to avoid this problem [Tof . 

A more severe problem of the SHDM scenario is the 
spatial anisotropy of the expected signal predicted due to 
the non-central position of the Sun in our Galaxy In 
a recent paper 70] it was shown that the non-observation 
of anisotropy in the data of the SUGAR experiment ex- 
cludes any SHDM scenario at the 5cr level assuming that 
all events above 4 x 10 19 eV are from SHDM sources. For 
the extreme case where SHDM is responsible only for 
UHECRs above 6 x 10 19 eV, the annihilation scenario is 
disfavored at least at 99% CL by the SUGAR data, while 
decaying SHDM still has a probability of ~ 10% to be 
consistent with the SUGAR data. 

The SHDM scenario is therefore disfavored by present 
experimental data and will be finally tested by the Pierre 
Auger experiment in the near future. We will therefore 
here focus on topological defect models with p = 1. 
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FIG. 6: Neutrino fluxes per flavor predicted by the three nor- 
malizations of the TD model of Fig. Qj] compared to future 
experimental sensitivities. The line key is as in Fig. El 

Fig. shows the results for nix = 2 x 10 13 GeV, with 
B = 10~ 12 G, and the moderately low theoretical esti- 
mate from Ref. |34| . These parameters lead to optimistic 
neutrino fluxes for the maximal normalization consistent 
with all data. For detailed earlier discussions of extra- 
galactic top-down fluxes see Refs. [72TF^| . 

Fig. shows that already the improved upper limit on 
the true diffuse photon background by EGRET implies 
too small a UHECR flux compared to the AGASA ex- 
cess at energies E > 10 20 eV. The parameters used in 
the figure represent the "best fit point" for this model; 



in particular for all other masses mx the disagreement is 
more severe. The new EGRET upper limit thus strongly 
disfavors extragalactic top-down scenarios. In addition, 
independently of this problem of overproduction of GeV 
7— rays, a non-observation of TD neutrinos by 2006 will 
rule out the possibility that extragalactic top-down mech- 
anisms significantly contribute to the UHECR flux, as 
can be seen from comparing Fig. with Fig.[S] 



V. NEUTRINO FLUX IN SCENARIOS WITH 
NEW HADRONS AS UHECR PRIMARIES 
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FIG. 7: Neutrino fluxes per flavor predicted by scenarios 
where UHECRs are explained as new hadrons produced as 
secondaries of accelerated protons with m = 3, z max = 2, 
compared to expected experimental sensitivities. The solid 
line is for a flux of primary protons peaked at E = 10 21 eV 
and the dashed line is for a primary proton flux oc E~ 2 up to 
-B max = 10 22 eV. The line key is as in Fig.0 

Supersymmetric (SUSY) models with a strongly inter- 
acting particle as lightest supersymmetric particle (LSP) 
or next-to- lightest SUSY particle (NLSP) are very inter- 
esting for explaining the AGASA excess above the GZK 
cutoff. Hadrons containing a gluino were first suggested 
by Farrar as UHECR primaries 0> E3 • This model with 
a light gluino together with a light photino as cold dark 
matter candidate is meanwhile excluded |76l l77| . How- 
ever, more general models with a light gluino or a light 
sbottom quark are still viable. 

In a recent paper |7Sj a model- independent, purely 
phenomenological approach was developed. Since the ob- 
served extensive air showers (EAS) are consistent with 
simulated EAS initiated by protons, any new primary 
proposed to solve the GZK puzzle has to produce EAS 
similar to those of protons. Experimentally still allowed 
are photons as UHECR primaries: at 90% CL, ~ 30% of 
UHECRs above ~ 10 19 eV can be photons However, 
the simplest possibility consistent with EAS observations 
is to require that the new primary is strongly interacting. 
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The requirements of efficient production in astrophysi- 
cal accelerators as well as proton-like EAS in the atmo- 
sphere ask for a light hadron, < 3 GeV, while shifting the 
GZK cutoff to higher energies results in a lower bound 
for its mass, > 1.5 GeV [7!j. From these requirements 
general conditions on the interactions of new UHE pri- 
maries were derived. The production of new hadrons in 
astrophysical objects was investigated in Ref. |78|. It was 
found that proton-proton collisions in astrophysical ac- 
celerators cannot produce sufficiently high fluxes of new 
primaries without contradicting existing measurements 
of photon |23J and neutrino fluxes [3j,|4j,|5j. In contrast, 
for a light shadron with mass < 3 GeV and the astro- 
physically more realistic case of UHE proton collisions 
on optical/infrared background photons there is no con- 
tradiction with existing limits. Also, the required initial 
proton energy is not too extreme, E < 10 21 eV, which 
may be achieved by astrophysical acceleration mecha- 
nisms. The only essential condition for the sources is 
that they should be optically thick for protons in order 
to produce these new hadrons. This condition applies to 
all models with new particles produced by protons. 

One of the important features of scenarios with new 
hadrons, and of any model in which the production cross 
section tr P7 ^s of a new particle S is much smaller than 
the total proton-photon cross section a pll is the high 
flux predicted for secondary high-energy neutrinos. This 
neutrino flux is approximately i r CR <T p7 / Cp-y^s m terms 
of the maximal contribution of S particles to the observed 
cosmic ray flux, Fen — (E/10 20 eV) -2 eV/(cm 2 s sr). 

Fig. [7] shows that for a primary proton flux oc -E~ 2 , 
this model can be restricted already with 3 years of 
AMANDA-II data. A non-observation of neutrinos by 
2006 will make it impossible to render the production 
cross section of the new hadrons consistent with existing 
limits in these scenarios. 



VI. THE Z-BURST SCENARIO 

In the Z-burst scenario UHECRs are produced by Z- 
bosons decaying within the distance relevant for the GZK 
effect. These Z-bosons are in turn produced by UHE neu- 
trinos interacting with the relic neutrino background |80| . 
If the relic neutrinos have a mass m u , Z-bosons can 
be resonantly produced by UHE neutrinos of energy 
E v ~ Mf /(2ra„) ~ 4.2 x f 21 eV (eV/m„). The required 
neutrino beams could be produced as secondaries of pro- 
tons accelerated in high-redshift sources. The fluxes pre- 
dicted in these scenarios have recently been discussed in 
detail in Refs. H^,|^|. In Fig. |8jwe show an optimistic 
example taken from Ref. ^ij^ . As in Refs. [32l l8lj no lo- 
cal neutrino over-density was assumed. The sources are 
assumed to not emit any 7— rays, otherwise the Z-burst 
model with acceleration sources overproduces the diffuse 
GeV 7— ray background [^J. We note that no known 
astrophysical accelerator exists that meets the require- 
ments of the Z-burst model [3^, ■ 
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FIG. 8: Flux predictions for a Z-burst model averaged over 
flavors and characterized by the injection parameters z m i n = 
0, 2 m ax = 3, a = 1, m = 0, fi max = 3 x f0 22 eV in Eq. Q for 
neutrino primaries. The sources are assumed to be exclusive 
neutrino emitters. All neutrino masses were assumed equal 
with m„ = 0.33 eV and we again assumed maximal mixing 
between all flavors. The line key is as in Fig. 



However, a combination of new constraints allows 
to rule out the Z-burst mechanism even for pure neu- 
trino emitting sources: A combination of cosmological 
data including the WMAP experiment limit the sum of 
the masses of active neutrinos to < I eV [83| . Solar 
and atmospheric neutrino oscillations indicate that in- 
dividual neutrino masses are nearly degenerate on this 
scale [3^] , and thus the neutrino mass per flavor must sat- 
isfy m v < 0.33 eV. However, for such masses phase space 
constraints limit the possible over-density of neutrinos in 
our Local Group of galaxies to < 10 on a length scale 
of ~ I Mpc [ijl . Since this is considerably smaller than 
the relevant UHECR loss lengths, neutrino clustering will 
not significantly reduce the necessary UHE neutrino flux 
compared to the case of no clustering. For the maxi- 
mal possible value of the neutrino mass m v ~ 0.33 eV 
the neutrino flux required for the Z-burst model is only 
in marginal conflict with the FORTE upper limit as 
shown in Fig. [SJ For all other cases the conflict is consid- 
erably more severe. Also note that this argument does 
not depend on the shape of the low energy tail of the pri- 
mary neutrino spectrum which could thus be even mo- 
noenergetic, as could occur in exlusive tree level decays of 
superheavy particles into neutrinos |8o| . However, in ad- 
dition this possibility has been ruled out by overproduc- 
tion of GeV 7— rays due to loop effects in these particle 
decays |86j . 

As was discussed in Ref. [3^, the Z burst scenario in- 
volving normal astrophysical sources producing neutrinos 
and photons by pion production within the source were 
already ruled out by the former EGRET limit. 
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VII. NEUTRINO FLUXES IN SCENARIOS 
INVOLVING ACTIVE GALACTIC NUCLEI 
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FIG. 9: Neutrino flux per flavor predicted for the AGN model 
from Ref . [87l ] for a uniform distribution of blazars (no redshift 
evolution) . The position of the peak is governed by the initial 
proton distribution. The normalization is determined by the 
amplitude of the accompanying 7— ray flux to the new diffuse 
EGRET flux estimate. The line key is as in Fig. 

Active Galactic Nuclei (AGN) are very promising sites 
of particle acceleration. Though there is no direct evi- 
dence of proton acceleration in these objects, according 
to the Hillas condition -E max ~ qBR for the maximal 
energy, where q is the charge, B is the magnetic field 
and R is the linear size of the acceleration region, AGN 
cores, jets or hot spots can be sites for acceleration of 
UHECRs up to the highest energies E ~ 10 19 - 10 21 eV. 
Once accelerated, protons can escape the AGN freely or 
can lose part of their energy in interactions with back- 
ground protons and photons. High energy neutrinos can 
be produced in AGN via pion production by accelerated 
protons. 

The total power of a given object in neutrinos can be 
related to the total power in 7— rays of MeV-GeV en- 
ergies. However, such a connection is not straightfor- 
ward because the observed high energy 7— rays can be 
produced by several mechanisms not involving neutrino 
production, for example proton synchrotron radiation in 
magnetic fields and inverse Compton scattering of MeV- 
GeV electrons. Moreover, the spectrum of neutrinos from 
AGN are even more difficult to predict than the spectrum 
of cosmogenic neutrinos. The reason is that besides the 
unknown model-dependent spectrum of primary protons, 
the spectrum of background photons is also not known 
and in general is model dependent. Also, AGN are di- 
vided into subclasses with different properties of the ob- 
served photon spectrum. Many of those properties do not 
directly relate to the possible neutrino spectrum. This 
means that it is very difficult to predict the space dis- 
tribution of those AGN which contribute to the neutrino 



flux from the distribution of AGN subclasses. At least 
these distributions can be very model dependent. 

Due to the above complications we suggest here a phe- 
nomenological approach to the prediction of neutrino 
fluxes. Within this approach the neutrino flux in most 
of AGN models can be approximately characterized by 
three parameters, namely the amplitude, width and po- 
sition of the peak (plateau) in the differential spectrum. 
The position of the peak (plateau) is related to the spec- 
trum of background photons. The combination of ampli- 
tude and width defines the total power in neutrinos which 
can be related to the total power in MeV-GeV photons 
produced in the same pion production reactions. Exper- 
imental bounds on neutrino fluxes can be converted to 
constraints on model parameters. 

A similar approach can be used for most existing mod- 
els which predict neutrino fluxes from AGN, see for exam- 
ple Ref. 123 ■ The only difference would be the connection 
of the phenomenological parameters of the observed neu- 
trino spectrum to the physical AGN parameters in the 
given model. 

As an example, we will use the model of 7— ray pow- 
ered jets of Ref. [87]]. In this model the high energy 7- 
rays are produced by accelerated protons interacting with 
the ambient photon fields (supplied, for example, by the 
accretion disk around the massive black hole) through 
photo-meson processes. At the same time those protons 
produce neutrinos which are emitted in the direction of 
the jet. Therefore, this model predicts a high neutrino 
flux comparable in power with the 7-ray flux. The de- 
tailed numerical simulations of proton acceleration in the 
central engine of the AGN [8^ show that the collimated 
jet of almost mono-energetic high energy protons (lin- 
ear accelerator) can be created in the electro-magnetic 
field around the black hole and the energy of those pro- 
tons can be converted into photons and neutrinos, while 
protons can be captured inside the source. The nucleon 
flux leaving the AGN is well below the observed cosmic 
ray flux in this scenario. Furthermore, since all nucleons 
leaving the source are well below the GZK cutoff, there 
is no cosmogenic contribution to the neutrino flux from 
these sources. 

Fig. [5] shows a typical prediction for the diffuse neu- 
trino flux in this model. The neutrino flux is maximized 
in such a way that the accompanying photon flux satu- 
rates the new EGRET bound. As seen from Fig. al- 
ready now AMANDA-H and RICE data start to restrict 
AGN models. Three years of AMANDA-II data will sig- 
nificantly restrict the parameter space of those models. 
This will also restrict the contribution of ir° production in 
AGN to the extragalactic diffuse 7— ray flux at EGRET 
energies. 

In the AGN model discussed above, blazars would be 
seen by neutrino telescopes as point-like sources with 
neutrino fluxes which are smaller or of the same order as 
the photon flux emitted by these same sources and which 
are detectable by 7-ray telescopes. The most probable 
sources were discussed in Ref. l90l. 
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VIII. CONCLUSIONS 

Based on our transport code we reconsidered neu- 
trino flux predictions and especially their maxima con- 
sistent with all current data on cosmic rays and updated 
upper limits on neutrino fluxes and the diffuse extra- 
galactic GeV 7— ray background. We discussed predic- 
tions for fluxes of cosmogenic neutrinos produced through 
pion production of UHECRs during propagation, and for 
fluxes produced by AGN. We showed that extragalac- 
tic top-down scenarios can not contribute significantly to 
the observed ultra-high energy cosmic ray flux if standard 
evolution histories and injection spectra are used. The Z- 
burst mechanism where ultra-high energy neutrinos pro- 
duce cosmic rays via interactions with relic neutrinos is 
ruled out except if cosmological neutrino mass limits are 
invalid and/or if relic neutrinos cluster more strongly 
than expected based on standard phase space principles. 
Only the case of maximal neutrino mass m„ ~ 0.33 eV 
consistent with large scale structure and CMB observa- 
tions is only moderately excluded. 

The fact that a good part of the speculative scenarios 



of UHECR origin are now ruled out makes the enigma of 
the highest energy particles an even more exciting subject 
of study in our opinion. 
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